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ABSTRACT 



We have investigated the variability of the binary X-ray pulsar, SMC X-l, in data from 
several X-ray observatories. We confirm the ~ 60-day cyclic variation of the X-ray flux in the 
long-term monitoring data from the RXTE and CGRO observatories. X-ray light curves 
and spectra from the ROSAT, Ginga, and ASCA observatories show that the uneclipsed 
flux varies by as much as a factor of twenty between a high-flux state when 0.71 second 
pulses are present and a low-flux state when pulses are absent. In contrast, during eclipses 
when the X-rays consist of radiation scattered from circumsource matter, the fluxes and 
spectra in the high and low states are approximately the same. These observations prove 
that the low state of SMC X-l is not caused by a reduction in the intrinsic luminosity of the 
source, or a spectral redistribution thereof, but rather by a quasi-periodic blockage of the 
line of sight, most likely by a precessing tilted accretion disk. In each of two observations 
in the midst of low states a brief increase in the X-ray flux and reappearance of 0.71 second 
pulses occurred near orbital phase 0.2. A similar brief flux increase near orbital phase 0.2 
was observed during a low-state observation that did not have sufficient time resolution to 
detect 0.71 second pulses. These brief increases result from an opening of the line of sight 
to the pulsar that may be caused by wobble in the precessing accretion disk. The records 
of spin up of the neutron star and decay of the binary orbit are extended during 1991-1996 
by pulse-timing analysis of ROSAT, ASCA, and RXTE PCA data. The pulse profiles in 
various energy ranges from 0.1 to > 21 keV are well represented as a combination of a pencil 
beam and a fan beam. Finally, there is a marked difference between the power spectra of 
random fluctuations in the high-state data from the RXTE PCA below and above 3.4 keV. 
In the higher energy range the spectrum has a sharp break at 3.3 Hz, with fitted power 
law indices of 0.45 and 1.76 below and above the break. No break is evident in the power 
spectrum below 3.4 keV, and the fitted power law index is 0.51. In both spectra there is a 
positive deviation from the fitted power law around 0.06 Hz that may be QPO. 



1. INTRODUCTION 

SMC X-l, located in the Small Magellanic Cloud (SMC), is the most luminous and rapidly spinning 
of the persistent accretion-powered X-ray pulsars. Its X-ray flux pulses with a 0.71 second period, it 
is eclipsed by its B0 supergiant companion for ~ 15 hours of the 3.892-day orbit, and its uneclipsed 
X-ray flux exhibits aperiodic variations on time scales from tens of milliseconds to months. Its rotation 
rate is increasing on a time scale of 10 3 yr, and its orbit is decaying on a time scale of 10 6 yr. With 
a well-determined distance, small amount of intervening matter, and extreme values of luminosity and 
spin rate, SMC X-l is a well-placed and particularly interesting object for detailed studies. 

X-rays from "an extended region or set of sources" in the SMC were discovered in a rocket observation 
by Price et al. (1971). Leong et al. (1971) identified a discrete source in the "wing" of the SMC in Uhuru 
observations and named it SMC X-l. Further analysis of these and additional Uhuru observations 
by Schreier et al. (1972) revealed the eclipse phenomena which established its binary nature. The 
identification of its optical counterpart as the B0I star Sk 160, suggested by Webster et al. (1972), 
was confirmed by Liller (1973) who observed an optical intensity variation with the same period as 
that of the X-ray eclipses. Lucke (1976) discovered the X-ray pulsations in rocket and Apollo satellite 
observations. Accurate values of the orbital elements of SMC X-l and rates of spin- up due to accretion 
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torque were derived by pulse-timing analysis of SAS 3 data by Primini, Rappaport, & Joss (1977) and of 
Ariel V data by Davison (1977). Measurements of the velocity curve of Sk 160 by optical spectroscopy 
( [Hutchings et al. 1977| ; [Reynolds et al. 1993| ), combined with the velocity curve of the pulsar, yielded 
narrow limits on the masses of both components of the binary system. An X-ray burst lasting ~ 80 s 
and aperiodic variability on time scales ranging from tens of milliseconds to hours were reported from 
EXOSAT observations by Angelini, Stella, & White (1991). Levine et al. (1993) (hereinafter LRD) 
determined the rate of orbital decay in a comprehensive analysis of the eclipse center times derived 
by pulse-timing analysis of data from various satellite observations through 1989, including their own 
observations with the Ginga satellite. Their publication includes a discussion of the implications of the 
observed decay for the orbital dynamics of the system. 

The circumstellar environment of the Sk 160/SMC X-l system has been the subject of several inves- 
tigations. Van Paradijs & Zuiderwijk (1977) attributed anomalies in five-color photometry of Sk 160 
to the optical emission from an accretion disk around the neutron star, and Tjemkes, Zuiderwijk, & 
van Paradijs (1986) interpreted the optical light curve in terms of ellipsoidal variations, disk emission, 
and X-ray heating effects. Hammerschlag-Hensberge, Kallman, & Howarth (1984) found evidence in 
UV spectra obtained from IUE of the influence of X-ray illumination on the stellar wind of Sk 160 
in the form of a deviation from spherical symmetry of the outflow and an anomalously low terminal 
velocity compared to stars of similar spectral type. Blondin & Woo (1995) modeled the wind dynamics 
of the Sk 160/SMC X-l system in three-dimensional hydrodynamical calculations which took account 
of the effects of X-ray ionization. Model light curves and spectra derived by Monte Carlo calculation of 
X-ray propagation through the stellar wind conformed well to the light curves and spectra of SMC X-l 
derived from Ginga and ROSAT observations, and in particular with the data obtained during eclipse 
transitions when the X-ray line of sight passes close to the surface of Sk 160 (Woo et al. 1995). 

Extended periods of very low X-ray flux of the uneclipsed source have been reported in observations 
with Uhuru ( [Schreier et al. 19?2j ), Copernicus ( [luohy fc Rapley 1975| ), Ariel V ( |Cooke 1978| ), and COS- 
ES (|Bonnet-Bidaud fc van der Klis 1981| ). An extended period of exceptionally high flux, beginning in 
1970 September and lasting ~ 100 days, was reported from a three-year record of observations by the 
Vela 5B monitor ( Whitlock fc Lochner 1994 ). Gruber & Rothschild (1984), in particular, reported a 
large amplitude high-low state variation in the X-ray flux in data from three ~ 80 day observations with 
the UCSD/MIT instrument on HE AO 1 in 1977 and 1978. They noted that the observed variations 
could be ascribed to either a quasi-periodic cycle with a duration of about sixty days or to a band- 
limited red noise process of unspecified origin; a strictly periodic variation was not consistent with their 
observations. They also suggested two possible physical mechanisms for a cyclic variation — intrinsic 
luminosity variation due to unstable mass transfer or periodic occultation by a precessing tilted accretion 
disk. 

In this paper we address several topics concerning the variability of SMC X-l and the evolution 
of the Sk 160/SMC X-l system. The observations and data reduction are described in §2. In §3, 
we present long-term monitoring data from the All-Sky Monitor (ASM) on the Rossi X-ray Timing 
Explorer (RXTE) and the Burst and Transient Source Experiment (BATSE) on the Compton Gamma- 
ray Observatory (CGRO) that confirm the ~ 60-day cyclic variability of the X-ray flux (brief reports of 
this result were made by Zhang et al. (1996) and Levine et al. (1996)). We then present an analysis of the 
X-ray light curves and spectra derived from various pointed-mode observations by the ASCA, ROSAT, 
Ginga, and RXTE observatories for the purpose of determining the cause of the cyclic variability. In 
§4, the records of spin-up of the neutron star and decay of the binary orbit are extended from mid 1991 
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to early 1996 with results derived by pulse-timing analysis of data from ROSAT, ASCA and RXTE. 
Pulse profiles in several energy ranges from 0.1 to > 21 keV are derived from these data, and modeled 
by a pulsar radiation pattern consisting of a pencil beam and a fan beam. Finally, in §5 we compare 
the power density spectra of the random fluctuations of the high-state flux in the high and low-energy 
portions of the data recorded by RXTE. 



2. OBSERVATIONS AND INITIAL DATA REDUCTION 

Bright X-ray sources have been monitored with the BATSE instrument in the 20-100 keV range since 
April of 1991 and with the ASM instrument in the 1.5-12 keV range since 20 February 1996. The 
BATSE instrument has been described by Fishman et al. (1994), and the ASM instrument by Levine 
et al. (1996). Source intensities are derived from BATSE data by fitting the predicted changes in count 
rates caused by earth occultations to the measured rates ( |Harmon et al. 1992|) . Intensities of known 
sources are derived from ASM data by fitting predicted shadow patterns cast by coded masks to observed 
patterns recorded by position-sensitive proportional counters. 

The Ginga data were obtained with the Large Area Counter (LAC) comprised of 8 proportional gas 
counters with a peak total effective area of 4000 cm 2 , a mechanically-collimated field of view of 1.1° x 
2.0° FWHM, an energy range from 1.5 to 36 keV, and an energy resolution of 20% FWHM at 5.9 keV. 
The Ginga observatory has been described by Makino (1987) and details of the LAC by Turner et 
al. (1989). 

Data from four ROSAT imaging observations were obtained with the Position Sensitive Propor- 
tional Counter (PSPC) in the focal plane of the X-ray telescope (XRT). The ROSAT mission has been 
described by Triimper (1983) and the XRT-PSPC system by Pfeffermann et al. (1987) and Aschen- 
bach (1988). The XRT-PSPC had an angular resolution of ~ 25", and was sensitive in the energy range 
from 0.1 to 2.5 keV. We denote the four ROSAT observations with the labels "1", "2", "3", and "C". In 
the numbered observations, SMC X-l was close to the center of the imaged field of view. Observation 
"C" was part of a survey of the SMC ( [Kahabka fc Pietsch 1995 ); in that observation SMC X-l was 
~40' off axis. 

The ASCA imaging data were obtained with the Solid-state Imaging Spectrometer (SIS) and the Gas 
Imaging Spectrometer (GIS) which are sensitive in the energy range from approximately 0.3 to 12 keV. 
The ASCA observatory has been described by Tanaka, Inoue, & Holt (1994). 

The pointed-mode RXTE data were obtained during the in-orbit checkout (IOC) period of operation 
with the Proportional Counter Array (PCA) which has a peak effective area of 7000 cm 2 , a mechanically- 
collimated 1° FWHM field of view, 18% energy resolution at 6 keV, and an energy range from 2 to 60 
keV QJahoda et al. 1996Q . 



Table |IJ summarizes the times and exposures of the pointed observations and the "state" of SMC X-l 
in the high-low cycle as defined below. 



2.1. Background Subtraction of Pointed Observations 

Background rates in the source region of the images produced by the ROSAT and ASCA observatories 
were estimated from the count rates in separate nearby regions. For ROSAT observations 1, 2, and 3 
we chose the source region to be a circle of radius 150" and the background region to be an annulus 
of inner radius 150" and outer radius 300". For the poorly focused ROSAT C observation, we chose 
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the source region to be a circle of radius 262.5" and the background region to be an annulus with radii 
of 262.5" and 487.5". We excluded from both the background and source regions a circular region of 
radius 150"(188" for the off-axis pointing) centered on the foreground coronal star HD 8191 - a weak 
soft X-ray source. Also, we excluded a circular region of radius 150" around the transient source RX 
J0117.6-7330 (Clark et al. 1997) in the two observations (2 & 3) in which it appeared. 

In the GIS image of the ASCA observations, we chose the source region to be a circle of radius 
383" and the background region to be an annulus of radii 383" and 766". In observation 2, a circular 
region of radius 150" around HD 8191was removed from both the source and background regions. In 
the SIS image, we chose the source region to be a square 380" on a side centered on SMC X-l and used 
the remainder of the chip as the background region. HD 8191 was not in the field of view of the SIS in 
these observations. 

Background corrections to the Ginga data were made according to the algorithm of Hayashida et 
al. (1989). Background corrections to the RXTE data were derived from the earth-occulted count rates 
which are listed in Table |2] together with the on-target rates in seven energy ranges. We note that 
the ROSAT images show several faint and soft sources near SMC X-l that were within the Ginga and 
RXTE fields of view. No corrections to the background rates were made for these sources because their 
total contribution to the Ginga count rate was less than a few percent of the SMC X-l eclipsed count 
rate ( |Woo et al. 1995| ), and less than one percent of the RXTE rates. 

3. QUASI-PERIODIC HIGH-LOW FLUX VARIATION 
3.1. Long-Term Monitoring Observations 

Figure 1 displays the X-ray light curves of SMC X-l measured by the ASM and BATSE instruments 
during times when the pulsar was uneclipsed and in the orbital phase range between 0.1 and 0.9 relative 
to eclipse center. The ASM data are binned in intervals of the 3.89-day orbital period, and the BATSE 
data in five-day intervals. 

A cyclic flux variation on a time scale of ~ 60 days is clearly evident in the ASM light curve. To 
quantify the duration of each cycle, which appears to be variable, we defined times of the beginning and 
end of each high state to be when the count rate increased/decreased to 2 counts s _1 . We determined 
these times by interpolating between successive points with values below and above, or above and below, 
2 counts s _1 . The times between successive rises and falls defined in this way are plotted against cycle 
number in Figure 2 where they are seen to range from below 50 to above 60 days with a generally 
decreasing trend. The ASM light curve shows that the flux often increased relatively quickly at the 
beginning of each high state and decreased more slowly at the end. 

To test for the presence of cyclic variability in the BATSE data, we divided the data into four subsets 
of equal duration and folded each subset into 5 phase bins with trial periods between 30 and 100 days. 
In this process, the mean flux and its standard deviation in each phase bin were calculated using the 
inverses of the variances of the individual measurements as weights. At each trial period, we calculated 
the reduced x 2 statistic of the light curve assuming that no signal was present. The values of the 
reduced \ 2 statistic are plotted in Figure 3 for each quarter of the BATSE light curve and for the entire 
BATSE light curve. Also shown are the results of the same analysis applied to the entire ASM light 
curve. Leahy et al. (1983) showed that if a sinusoidal signal with a fixed period, P, is present in the 
light curve, such a folding analysis will show a peak centered at that period with a FWHM of P 2 /T, 
where T is the length of the observation. In Figure 3 the horizontal bars show the expected FWHM 
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of a sinusoidal signal with a 55-day period. We note that the positions of the peaks shift and some of 
peaks are substantially wider than the horizontal bars. These features, together with the drift of the 
time since last increase/decrease in the ASM light curve (Figure 2), are strong evidence for the absence 
of a long-term coherent periodicity in these data sets. Instead, the observed cyclic variation can be 
best characterized as quasi-periodic with recurrence times that wander between ~50 and ~60 days. 
Attempts to examine the variation of the cycle more closely in the BATSE light curve by smoothing 
and by correlation with a sine wave were inconclusive, presumably due to the low signal-to-noise ratio 
in the data. 

3.2. Short-Term Pointed Observations 

Background-subtracted count rates, derived from the Ginga, ROSAT, ASCA, and RXTE PCA ob- 
servations, are plotted in Figure 4 against orbital phase. We multiplied the background-subtracted 
count rates of the off- axis ROSAT observation by a factor of 1.7 to compensate for the reduction in 
effective area relative to the on-axis observation. So that the Ginga and ASCA count rates may be 
directly compared, we show the count rates from the pulse height channels which correspond to the 
energy range 1-10 keV. For the ASCA GIS, these are PI channels 85 through 850. For the Ginga LAC 
these are channels 3 through 18 in the MPC-2 mode data and 1 to 4 in the MPC-3 mode data. 

With the exception of the Ginga observation, which occurred before the beginning of the BATSE light 
curve, all of the pointed observations discussed in this paper occurred during the time of the BATSE 
observations, and before the RXTE ASM observations. Unfortunately, the signal-to-noise ratio of the 
BATSE light curve is not large enough to allow certain determination of the phases in the flux cycle at 
which the individual pointed observations occurred. Furthermore, because of the quasi-periodic nature 
of the cycle, it is not possible to extrapolate backwards from the ASM light curve. However, the ASM 
light curve shows that the X-ray flux averaged over the uneclipsed portion of each binary orbit rises and 
falls smoothly over each cycle. Therefore, we are confident that the observed out-of-eclipse flux of the 
source, as it is revealed in the light curves, is a reliable indicator of the state of the source. 

The contrast between the high and low states of SMC X-l is most clearly demonstrated by the 
various ROSAT observations. Where observation 1 and the region C observation overlap in orbital 
phase outside of eclipse, they differ in (corrected) count rate by a factor of at least 20. We therefore 
identify ROSAT observation 1 as having occurred during the high state and the region C observation as 
having occurred during the low state. The light curve of ROSAT observation 2 behaves like that of the 
region C observation where the two overlap so we identify this observation as having occurred during 
the low state as well. Outside of eclipse, ROSAT observation 3 shows a count rate comparable to that 
of observation 1, so we identify it as having occurred during the high state. 

The two ASCA observations also differ by a factor of about twenty in their out-of-eclipse count rates, 
so we assume ASCA observation 1 occurred during the high state and ASCA observation 2 occurred 
during the low state. However, the decline in the count rate in observation 1 beginning near orbital 
phase 0.5 may be a transition to the low state. In the Ginga observation, the out-of-eclipse behavior 
and count rate is comparable to ROSAT observation 1 and ASCA observation 1, so we conclude it also 
occurred during the high state. The brief RXTE PCA observation near orbit phase 0.5 clearly occurred 
during a high state. 

Interesting features of the low-state light curves derived from the ROSAT observations 2 and C are 
the increases in the count rates near orbital phase 0.25. In both cases, 0.71 second pulsations are present 
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during the increase. A similar turn-on may be seen in the light curve (not shown herein) derived from 
the ROSAT All-Sky Survey ( [Kahabka fc Pietsch 19961) , although those data do not have sufficient time 
resolution to permit a test for the presence of pulsations. ROSAT observation 2 ends near phase 0.3 with 
this increase. The ROSAT C observation has a gap from phase 0.3, where the count rate is increasing, 
to phase 0.58 (plotted at phase -0.42), where the count rate is lower than before. 

In Figure 5 we plot the photon number spectra of SMC X-l from all the observations during the 
in-eclipse portions of the binary orbit (orbital phase -0.07 to +0.07). For the purpose of comparing 
these spectra, we fitted all of them simultaneously with the following model spectrum that consists of 
the sum of a black-body spectrum and a power law with a high energy cutoff together with interstellar 
absorption: 

1(E) = ex P [-a(E)N u }[f hh (E) + f pl (E)\, (1) 

where 

f hh (E) = I hh (E/E hh ) 2 (e - l)[exp(E/E hh ) - 
f pl (E) = I pl (E/lkeV)~ a x { { _ E _ Ec} 




I ex Pi— Et 

Jbb is the the intensity of the blackbody component at E = E^, I p \ is the intensity of the power law 
component at E = 1 keV, and cr(E) is the photoelectric absorption cross section of cold matter ( |Morrison 



& Mc Cammon 1983] ). The results of the fit are given in Table ^. The convolutions of this model with 
each of the instrument responses are plotted in Figure 5 along with the corresponding observed spectra. 
The instrument response for the ROSAT C observation includes the reduced efficiency at the off- axis 
location. 

The high-low state identifications we have inferred from the observed fluxes leave no doubt that the 
spectra displayed in Figure 5 include instances from both the high and low states during the eclipse when 
the only X-rays detected are those that have been scattered from circumsource matter. Since these five 
spectra show no more than a factor of two difference in flux at any X-ray energy, we conclude that the 
average illumination of the circumsource matter varied little between the high and low states. Moreover, 
all five spectra have similar shapes. Thus, the cause of the cyclic variation cannot be either a variation 
in the intrinsic luminosity of SMC X-l or a redistribution of the luminosity across the spectrum. We 
conclude that the cause of the cyclic variation between states of high and low flux is a quasi-periodic 
blockage of the line of sight. The same conclusion regarding the cause of the long-term cyclic variation 
in the X-ray flux of LMC X-4 has been drawn on the basis of similar evidence ( |Woo, Clark, fc Levine 
19950 



3.3. Motions of the Accretion Disk 



Shortly after discovery of the 35-day high-low state variation in Her X-l, Katz (1973) suggested that 
if the rim of the neutron star's accretion disc does not lie in the plane of the binary orbit, then the plane 
of the rim will precess due to the gravitational force of the companion star. Katz argued that if the rim 
is sufficiently tilted, it will periodically block the line of sight to the neutron star, causing the observed 
cyclic high-low variation. His expression for the forced precession period, Pf , of a ring of non-interacting 
particles can be written in the form 



Pi 



3 Pb U: 



3/2 



0.6g 2 / 3 + ln(l + g 1 / 3 ) 
0.49g 2 / 3 



3/2 



+ q 



2x1/2, 



cos j3)~ 



(2) 
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where Pb is the binary orbital period, Pl is the volume-averaged radius of the Roche lobe, P r is the ring 
radius, q is the ratio of the neutron star mass to the companion star mass, and j3 is the tilt angle of the 
ring with respect to the plane of the binary orbit. The quantity in square brackets is an approximate 
expression for the ratio of the binary separation to the Roche lobe radius (Eggelton 1983). 

Analyses of the optical light curves of Sk 160/SMC X-l indicate that the accretion disk radius is ~ 0.7 
to 1.0 times the Roche lobe radius QHowarth 1982| ; [Tjemkes, Zuiderwijk, fc van Paradijs 1986 ; |Khruzina 
& Ch frepashchuk 1987] ) . Setting R V /R L = 0.7, q = 0.0909 ( |Reynolds et al. 1993| ), P b = 3.892 days, and 
cos/5 ~ 1, we find the precession period predicted by equation || is 31 days or 0.56 times the typical 
observed cycle time of 55 days. The optical light curves of the two other X-ray binaries with clear quasi- 
periodic high-low state cycles, Her X-l and LMC X-4, also indicate that their accretion disks nearly fill 
the Roche lobes of their neutron stars ( |Gerend fc Boynton 1976 ; Qovaisky et al. 1984 ; Hccmskcrk fc van 



Parad pjs" 1989| ). The corresponding precession periods (assuming R v /Rl = 0.7 and cos (3 ~ 1) are 0.56 
and 0.35 times the observed high-low cycle periods of Her X-l and LMC X-4, respectively. Thus, in all 
three cases a simple model of an occulting rim of non-interacting particles predicts precession periods 
substantially less than the observed high-low cycle times. 

Numerous theoretical investigations of the effects of the accretion rate, viscosity, and X-ray illumina- 
tion on the creation of warp and tilt in the disk and on the precession rate have been carried out. Among 
these investigations we note in particular the work of Iping & Petterson (1990). They made numerical 
simulations with results that showed stable precessing disks with quasi-periods that could be adjusted 
to fit observations by appropriate choices of the accretion rate, the conversion efficiency with which 
accretion energy is converted into X-rays, and the viscosity parameters. Pringle (1996) confirmed their 
results with a more rigorous treatment, and showed that a warp instability due to radiation pressure will 
develop outside a critical radius which is ~ 3 x 10 8 cm for a neutron star. The radius of the Roche lobe 
in SMC X-l is ~ 2.9 x 10 11 cm so the accretion disk is clearly large enough to be unstable to warping. 
Thus it appears that a detailed model of a precessing tilted accretion disk could provide a quantitative 
explanation for the observed high- low cycle of SMC X-l. 

It is unlikely that random flares were the cause of the short-term flux increases observed near orbital 
phase 0.25 in the three low-states observations. The presence of 0.71 second pulses in the ROSAT 2 
and C observations shows that the increased flux came directly from the neutron star. It is also unlikely 
that the increase in the ROSAT C observation was the beginning of a high-state because a decrease 
followed soon after, and because the observation occurred only ten days after the source was observed 
to be in a high state in ROSAT observation 1. The ASM light curve shows that low states typically 
last for 15-20 days. Thus the ROSAT C observation probably occurred in the early part of a low state. 
These increases are reminiscent of the low-state turn-ons in Her X-l which occur near orbital phases 0.2 
and 0.7 (cf. Fig. 2 of Priedhorsky & Holt 1987), and which have been attributed ([Levine fc Jernig"an 



1982Q] to a wobble which may be a common feature of precessing accretion disks. Further observations 
of SMC X-l are needed to determine the regularity and nature of its low-state turn-ons. 



4. PULSE-TIMING ANALYSIS OF SPIN UP AND ORBITAL DECAY 

4.1. Observations and Analysis 

We carried out pulse timing analysis of the ROSAT, ASCA, and RXTE PCA data to extend previously 
published results on the spin up of the neutron star and decay of the binary orbit through 1996 January. 
The ROSAT and ASCA data consist of the energy, image position, and readout time of each photon 
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detected in an image of SMC X-l. The time resolution of the ROSAT data is a few ms, and the time 
resolution of the ASCA data is 1/16 s. The RXTE data consist of the counts accumulated by the PCA 
in rebinned time intervals of 1/64 s and five energy channels from 2 to > 21 keV. 

None of our six pointed observations span a sufficient portion of the SMC X-l orbit to permit an 
independent determination of all of the orbital parameters and the rate of change of the pulse frequency. 
The ROSAT 2 and C observations had such brief exposures to the pulsed flux that the data were only 
sufficient to constrain the pulse frequency when the orbital phase (measured from eclipse center for a 
circular orbit) was fixed at the value inferred from the other observations. We restricted the timing 
analysis of the other four observations to determinations of the intrinsic pulse frequency and orbital 
phase at the start of each observation, while the other orbital parameters were fixed at the values 
derived by LRD from observations of longer duration. From the orbital phases we derived the center 
times of the preceding eclipses. These were combined with the previous eclipse center times quoted by 
LRD in a redetermination of the orbital decay rate. 

To obtain initial estimates of the orbital phase, 4>o, and intrinsic pulse frequency, /, at the start of 
an observation, we searched a two-dimensional array of phase and frequency for the pair of values that 
maximized the fractional amplitude of the folded pulse light curve derived from the entire set of solar 
system barycentric X-ray arrival times corrected for orbital delay. The delay-corrected arrival time of 
the zth photon, t'^ was derived from its observed arrival time, £j, by iterative solution of the equation 

t'. 

t[ = ti — a x sin i cos(27r— + O ) (3) 

where the quantities a x sin i and a\ were fixed at the values found by LRD and listed in Table [5[ The 
change in these values, implied by the rate of orbital decay measured by LRD, had a negligible effect 
on the delays computed for the present data. 

We refined the estimated values of 0q and / by the following procedure aimed at obtaining a sequence 
of accurate pulse arrival time delays to which calculated orbital delays could be fitted by adjustment 
of the orbital phase. A portion of the delay- adjusted data was folded at the estimated frequency, and 
a 12-term Fourier series was fitted to the resulting pulse profile to form a template profile. The data 
were then divided into a number of subsets, and each subset was folded with the estimated frequency 
to obtain its pulse profile. Each profile was cross-correlated with the template to determine its relative 
phase shift due to errors in the estimated values of 0o an d /. From the phase shift of the jth profile 
we computed the difference between the actual corrected arrival time of a pulse near the mean arrival 
time, Tj, of the data subset and the expected corrected arrival time based on the estimated frequency. 
We call r(Tj) the sum of that difference and the orbital delay at Tj calculated with the estimated 0o- 
The quantity r{Tj) is equal to the sum of the true orbital delay and the time equivalent, (Tj — T )Sf/f, 
of the portion of the phase shift due to the error, 5f, of the estimated value of /. A least-squares 
fit of calculated values of r(Tj) to the set of measured values was made by adjustment of </> and Sf. 
This procedure was repeated, each time using improved estimates of the orbital phase and frequency 
to compute a new template and new values of r[Tj) until there was no further improvement in the 
least-squares fit, and no significant trend away from a zero shift in pulse phase between the template 
and the subset pulse profiles. 

In Figure 6 we show plots of the final pulse arrival time delays together with the fitted curves and 
deviations of the fits. The ROSAT 1 and 3 and the ASCA 1 observations yielded well-constrained 
values of the intrinsic pulse period P = 1/ f and eclipse center times. In spite of the high statistical 
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accuracy of the RXTE data, the shortness of the observation resulted in relatively large and correlated 
uncertainties in the eclipse center time and intrinsic pulse frequency. 

The quadratic function of orbit number defined by LRD was fitted to the set of eclipse center times 
which includes the four new ones derived in this work and the seven previous ones listed in their paper 
and reproduced in Table |j. The result is shown in Figure 7 as a plot of the measured deviations from 
a linear function of orbit number against time. The fitted values of the coefficients are listed in Table |5| 
together with the previous values of LRD with which they are entirely consistent. 

The pulse periods from this and previous studies are plotted against time in Figure 8 together with 
the deviations from a fitted quadratic function of time with the coefficients listed in Table ^|. We note 
that every measurement of the pulse period of SMC X-l has yielded a value lower than the previous 
measurement. Although the magnitude of the rate of change of the period appears to be decreasing, the 
absence of any observed episodes of torque reversal indicates that the accretion flow has never slowed 
enough to allow the inner edge of the accretion disk to retreat to the corotation radius where braking 
action would be expected to occur. 

4.2. Variation of the Pulse Profile with Time and Energy 

Figure 9 displays background-subtracted pulse profiles in five energy ranges from 0.1 keV to > 21 
keV derived from the ROSAT, ASCA, and RXTE PCA data. 

Striking features of the profiles are the simplicity and smoothness of the profiles, and the steady 
increase with energy in the relative amplitude of the narrower component of the double pulse. In 
marked contrast, the profiles of various other X-ray pulsars such as Vela X-l and LMC X-4 are complex 
and vary widely over the same energy range. Such complexity may be caused by a multi-pole structure 
of the magnetic field that guides the accretion flow to the surface of the neutron star, giving rise to 
a multiplicity of hot spots and complicated patterns of emission from multiple accretion columns. In 
SMC X-l it appears that the guiding field may be represented to a good approximation by a simple 
dipole inclined with respect to the spin axis. The different energy dependence of the two components, 
and the fact that the amplitude of the narrower one changes from less to more than that of the broader 
one precludes the possibility that they are caused by identical pencil beams emanating in opposite 
directions from the magnetic poles. Instead, the narrower pulse may be attributed to the pencil beam 
emanating from the accretion column above one pole, and the broader pulse to the fan beams emanating 
from the accretion columns above the two poles. 

A second prominent feature is the variability with time and energy of the narrower pulse. In the 
energy bands where data from ASCA in April 1993 and from RXTE in January 1996 can be compared, 
there is a substantial change in the relative amplitudes of the two pulses which appears to be due to 
a change in the flux in the narrower pulse. Similar variability of the narrower pulse was reported by 
LRD in Ginga data from 1987, 1988, and 1989. In each of the three data sets displayed in Figure 9 the 
amplitude of the narrower pulse increases with energy, starting near zero in the ROSAT data below 1 
keV, and rising above that of the wider pulse in the RXTE data above 21 keV. 

X-ray production and radiative transfer in magnetized plasma has been computed by Nagel (1981) for 
slab and cylindrical geometries. In one particular case of cylindrical geometry with the magnetic field 
parallel to the axis, which serves as a model of an accretion column, the resulting radiation pattern has 
the form of combined pencil and fan beams, i.e., the escaping ordinary photons form a pencil beam with 
a maximum near 0° with respect to the axis, and the extraordinary photons form a fan beam with a 
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maximum at 90°. It may be plausible, therefore to model the SMC X-l pattern as such a combination. 
If the broad pulse is identified with a fan beam, then its flat top can be attributed to the prolonged 
dwell near the minimum angle between the line of sight and the plane perpendicular to the dipole axis. 
If the narrow pulse is attributed to the pencil beam, its energy dependence and variability may be due 
to the energy dependence and variability of the opacity of the magnetized plasma through which the 
pencil beam propagates. 

In Figure 10 we display the observed pulse profile of background-subtracted data from 2 to 21 keV 
(histogram) together with a model profile (crosses) computed for a radiation pattern consisting of an 
isotropic component plus combined pencil and fan beams. The pencil beam flux is represented by a 
Gaussian function of the angle, p, between the line of sight and the dipole axis, the fan beam flux by a 
Gaussian function of 90° — p. We call G the inclination of the spin axis to the line of sight, and ^ the 
angle between the spin axis and magnetic dipole axis. Then 

p = cos -1 [sin Q cos(0 — O ) srn ^ + cos G cos \ff] . (4) 

where is the pulse phase. The predicted flux was computed by the equation 

m = Fr exp (- Jj) + F, exp + F„. (5) 

We set G to 60°, which implies equality between the the inclination of the spin axis and estimates of the 
orbit inclination, and ^ to 27°. This combination of angles achieves the condition of a dwell of the fan 
beam near the line of sight. The other parameters were adjusted by eye and are listed in Table ^ We 
note that the pencil beam in this model is similar in angular width to that in the theoretical model of 
Nagel (1981) mentioned above. On the other hand, the fitted fan beam is substantially narrower than 
the computed one. 



5. RANDOM INTENSITY FLUCTUATIONS 

Power spectra of the random fluctuations in the RXTE PCA data in two energy ranges are displayed 
in the right-hand panels of Figure 11. The spectra were computed as follows: the data were divided 
into 53 segments of 256 seconds duration, the power spectrum of each segment was computed with the 
normalization of Leahy et al. (1983), the Poisson noise component of 2 was subtracted, and the result 
was divided by the mean count rate of the segment to obtain fractional RMS normalized powers. The 
53 spectra were then averaged, and the powers due to five harmonics of the pulsed power, computed 
according to the method of Angelini et al. (1991), were subtracted. Finally, the resulting power spectra 
were logarithmically rebinned for display. The broken power laws and harmonic components that were 
fitted to the normalized spectra are displayed in the left-hand panels of Figure 11. The errors of 
the normalized powers before rebinning were all set equal to the expected RMS Poisson fluctuations 
computed as the quantity 

" M 

£(!/< 



2 

M 



1/2 

(6) 



where q is the count rate during the ith data segment, and M is the number of segments. 

A marked difference is apparent between the forms of the power spectrum in the low and high energy 
channels. The spectrum in the 2.0-3.4 keV range is fitted by a single power law with an index of 0.51; 
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the 3.4-21 keV spectrum conforms to a broken power law with indices of 0.45 and 1.76 with a sharp 
break at 3.3 Hz. Both spectra show a positive deviation from the power-law fit around 0.06 Hz which 
may be attributed to quasi-periodic oscillation (QPO). 

The break we find in the power spectrum above 3.4 keV is similar to those found by Takeshima (1994) 
in Ginga observations of several pulsars. He fitted the power spectra with broken power laws, and found 
breaks occurring near the fundamental pulsar frequencies. In his analysis of data from the second and 
third Ginga observations of SMC X-l listed in Table |], he found breaks at 2.93 Hz and 1.79 Hz. 

The difference between the power spectra that we find in the high and low-energy portions of the 
X-ray flux may be a clue to the effect of magnetic field intensity on the frequency spectrum of emission 
fluctuations. The flux in the low-energy channel may be dominated by radiation emitted from the 
accretion disk where the magnetic field is relatively weak. Most of the flux in the high-energy channels 
certainly arises from the accretion column and hot spot near the surface of the neutron star, where the 
magnetic field is strong. The strong field may suppress the fluctuations responsible for low-frequency 
red noise, thereby causing the abrupt decrease in the logarithmic slope of the power spectrum below 3.3 
Hz. 

6. Summary 

The principal results in this paper are as follows: 

The ASM light curve of SMC X-l confirms the presence of a high-low state cyclic variation of the 
X-ray flux, and shows that the cycle duration varied from ~ 50 to ~ 60 days. The BATSE data show 
that the cyclic variation is a persistent phenomenon. 

Results from pointed observations show that the uneclipsed flux of SMC X-l in the high state is more 
than an order of magnitude greater than in the low state. In contrast, the eclipsed flux, consists of 
X-rays scattered from circumsource matter, differs by no more than a factor of ~ 2 between the two 
states, and the spectrum shape does not change significantly. Thus, the high-low cyclic variation is not 
caused by either a change in the luminosity of the source, or a redistribution of the energy. Rather, it 
must be the result of a quasi-periodic blockage of the line of sight. 

In a simple model of occultation by a ring of non-interacting particles at the rim of a tilted accretion 
disk, the forced precession caused solely by the gravity of the companion star would have a period of 
about half the observed long-term cycle duration. Similar differences are found between the predicted 
precession periods for the simple model and the observed long-term cycle periods of Her X-l and LMC X- 
4. It has been shown that more detailed and realistic models of precessing accretion disks that take 
account of the effects of the accretion rate, accretion energy conversion efficiency, and viscosity in the 
disk, can be adjusted to cause model disks to tilt and precess with periods in agreement with the 
observed long-term cycles. We therefore conclude that the high-low state cycle in the X-ray flux from 
SMC X-l, like the similar cycles of Her X-l and LMC X-4, is the result of quasi-periodic occultation by 
a precessing tilted accretion disk. 

Timing analysis of the 0.71 second pulses shows a continuation of the decay of the binary orbit through 
1996 January at a rate consistent with previous results. 

The pulse period has continued its apparently monotonic decrease at a gradually decreasing rate. 

The pulse profile can be modeled with fair accuracy by a three component radiation pattern consisting 
of an isotropic component plus combined pencil and fan beams symmetric about an axis inclined at 
27° from the spin axis, with the spin axis inclined at 60° to the line of sight. 
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The power spectrum of the flux recorded by RXTE in energy channels below 3.4 keV is fitted by a 
single power law with index 0.51 over the frequency range from 0.01 to 32 Hz. In marked contrast, the 
power spectrum in energy channels above 3.4 keV is fitted by a broken power law with a sharp decrease 
in index from 1.76 above 3.3 Hz to 0.45 below. The low-energy flux may be dominated by radiation 
emitted from the accretion disk where the magnetic field is relatively weak. The high-energy flux is 
emitted primarily in regions of strong magnetic field which may suppress the emission fluctuations that 
cause the low-frequency red noise. 

In both the high and low-energy power spectra there is a positive deviation from the fitted power law 
around 0.06 Hz that may be QPO. 

We used the XSPEC X-ray spectral fitting program for our spectral analysis. Archival ROSAT and 
ASCA data was obtained with the use of the High Energy Astrophysics Science Archive Research Center 
Online Service, provided by the NASA/Goddard Space Flight Center. We thank Jane Turner and Ken 
Ebisawa of the user support teams for advice on reducing and analyzing ROSAT and ASCA data, and 
Edward Morgan of MIT for assistance in the use of RXTE data acquired during the in-orbit calibration 
(IOC) period of operation. We specially thank the referee, J. van Paradijs, for his careful reading and 
valuable suggestions for improvements in the manuscript. 
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TABLE 1 
Summary of Pointed Observations 



Instrument and ID 


Start Time-Stop Time 


Phases a 


Exposure (s) 


State 


ref 


Ginga 


1989 Jul 30.78-Aug 3.20 b 


0.34 - 0.21 


74240 


high 


1 


ROSAT 1 


1991 Oct 7.17-8.13 


0.48 - -0.29 


16634 


high 


2 


ROSAT C 


1991 Oct 16.31-19.98 


-0.18 - -0.24 


22173 


low 


3 


ROSAT 2 


1992 Sep 30.71-Oct 2.59 


-0.15 - 0.33 


8932 


low 


2 


ROSAT 3 


1993 Jun 2.99-4.20 


-0.13 - 0.18 


11864 


high 


2 


ASCA 1 


1993 Apr 26.94-27.97 


0.35 - -0.39 


34507 


high 




ASCA 2 


1995 Oct 18.79-19.88 


-0.17-0.11 


30628 


low 




RXTE 


1996 Jan 11.00-11.23 


0.47- -0.47 


14536 


high 





a Binary orbit phase range on the (-0.5,0.5) interval. All observations span less than one binary orbit. 
b SMC X-l was observed for more than one binary orbit. Only a subset of the data actually collected 
is described here. 

References.— (1) Woo, et al. 1995; (2) Clark, et al. 1997; (3) Kahabka & Pietsch 1996. 
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TABLE 2 

Average Count Rates in Seven Energy Ranges 
of RXTE Data 



Average Count Rates 





(counts s" 


- 1 ) 


Energy Range 


SMC X-l 


Earth 


(keV) 


+Background 


Occulted 


2.0-3.4 


43.5 


5.7 


3.4-5.8 


173 


5.0 


5.8-8.6 


173 


4.6 


8.6-21 


238 


21 


21-60 


74 


53 
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TABLE 3 

Fitted Values of the Spectral Function Parameters 



Parameter Value (la error) 



N n {10 20 cm" 2 ) 5.9 (4.2) 

E hh (keV) 0.25 (0.02) 

4b (10~ 3 photons cm- 2 s" 1 keV" 1 ) 9. (3) 

a -0.2 (0.1) 

J pl (10- 5 photons cm" 2 s" 1 keV" 1 ) 9.8 (3.1) 

E c (keV) 6.5 (0.4) 

E f (keV) 5.9 (0.4) 

Xl 3.96 
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TABLE 4 

Summary of Mid-Eclipse Times and Pulse Period Measurements 



U bservatory 


(MJrJJJ 


AT 
1 V 


r pulse [i-CT ) 
S 


Keierence 


Uhuru 


40963.99(2) 
41114.2 b 


-481 




Schreier et al. 1972 


Uhuru 




0.71748(26) 


Henry & Schreier 1977 


Aerobee 


41999.6 b 




0.7164(2) 


Yentis et al. 1977 


Copernicus 


42275.65(4) 
42613.8 b 


-144 




Tuohy & Rapley 1975 


Apollo-Soyuz 




0.7151(2) 


Yentis et al. 1977 


SAS 3 


42836.1828(2) 





0.71488585(4) 


Primini et al. 1977 


Ariel V 


42999.6567(16) 
43985.907 b 


42 


0.7147337(12) 


Davison 1977 


Einstein 




0.713684(32) 


Darbro et al. 1981 


COS B 


43116.4448(22) 


72 




Bonnet-Bidaud et al. 1981 


Ginga 


46942.47237(15) 


1055 


0.710592116(36) 


LRD 


Ginga 


47401.744476(7) 


1173 


0.7101406720(15) 


LRD 


Ginga 


47740.35906(3) 


1260 


0.709809901(32) 


LRD 


ROSAT 1 


48534.34786(35) 
48892.4191 b 


1464 


0.70911393(34) 


Present work 


ROSAT 2 


1556 


0.7088166(6) 


Present work 


ASCA 


49102.59109(82) 


1610 


0.7086390(9) 


Present work 


ROSAT 3 


49137.61911(50) 


1619 


0.7085992(5) 


Present work 


RXTE 


50091.170(63) 


1864 


0.70767331(8) 


Present work 



a Mid-eclipse time, equivalent to time when mean longitude I equals 7r/2 for a circular orbit (c.f. 
Deeter, Boynton, fc Pravdo 1981| ). 
b Set to predicted mid-eclipse times 
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TABLE 5 

Parameters of the Spin-Up and Orbital Decay 



Value (la uncertainty) 



Parameter Present Work LRD 



a x sini (lt-s) 




53.4876 


Eccentricity 




< 0.00004 


P pulse / P pulse (y r ) 


-4.52(2) x KT 4 a 


-5.30(2) x KT 4 


a b (MJED) 
ai b (days) 
02 b (days) 


42836.18278(20) 


42836.18277(20) 


3.89229090(43) 


3.89229118(48) 


-6.953(28) x 10~ 8 


-6.97(3) x 10~ 8 


Porb/Porb (J!" 1 ) 


-3.353(14) x 10~ 6 a 


-3.36(3) x 10~ 6 



Evaluated at MJED 50000. 

b Center time of the N th eclipse is given by t N = a + aiN + a 2 N 2 . 
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TABLE 6 

Coefficients of the Quadratic Fit to the 
Pulse Period Data 



Parameter Value 



Po 


0.707753 


s 


Pi 


-8.75 x 10- 6 


s day -1 


P2 


1.68 x 10~ n 


s day -2 


to 


2450000.5 


JD 


-Ppulsc 


= Po + Pi(t - t ) +P2(t 


-to) 2 - 



21 



Fig. 1.— CGRO BATSE and RXTE ASM (1.5-12 keV)light curves of SMC X-l. The time of the start 
of each of the pointed observations is indicated on this figure by R for ROSAT, A for ASCA, and X for 
RXTE. 

Fig. 2. — Time intervals between successive high/low transitions of SMC X-l. The transition times are 
obtained from crossings of 2 counts s" 1 in the ASM light curve. 

Fig. 3. — Deviations from the average count rates, given in terms of a reduced chi-squared statistic, 
of four sections of the BATSE light curve, of the entire BATSE light curve, and the ASM light curve 
when folded at various trial periods. The horizontal bars indicate the expected width (FWHM) of a 
peak corresponding to a sinusoidal variation with constant period. 

Fig. 4. — Count rates of SMC X-l as a function of orbital phase for the pointed observations. The 
ROSAT PSPC region C count rates are corrected for off-axis vignetting. Dashed lines at orbital phase 
±0.07 indicate the end of immersion and beginning of emersion of the eclipse transitions observed in 
the high state. 

Fig. 5. — Observed photon number spectra during eclipses (orbital phase in the range -0.07 to +0.07) 
and the computed spectra (solid lines) for a single fitted model. 

Fig. 6. — Pulse arrival-time delays plotted against orbital phase together with the fitted curves of 
predicted delays and residuals. 

Fig. 7. — Differences between the measured values of eclipse center times (7^/2 ) and values computed 
for a constant orbital period (P orb =3. 892188). The solid curve is the fitted function (see Table 5) which 
includes a linear decrease of the orbital period. 

Fig. 8. — Inferred intrinsic pulse period plotted against time together with a fitted quadratic and 
residuals. 

Fig. 9. — Pulse profiles in various energy bands derived from the ROSAT 1, ASCA 1, and RXTE data. 

Fig. 10. — Computed pulse profile (crosses) modeled by pencil and fan beams and fitted to the observed 
average 2-21 keV pulse profile (histogram) derived from the RXTE PCA data. 

Fig. 11. — Left panels: Fitted harmonic and continuum components of the power spectra in two energy 
bands of RXTE PCA data with times corrected for orbital motion. Right panels: Power spectra after 
subtraction of the harmonic components. 
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